Spontaneous radiation of a finite-size dipole emitter in hyperbolic media 
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We study the radiative decay rate and Purcell effect for a finite-size dipole emitter placed in a 
homogeneous uniaxial medium. We demonstrate that the radiative rate is strongly enhanced when 
the signs of the longitudinal and transverse dielectric constants of the medium are opposite, and 
the isofrequency contour has a hyperbolic shape. We reveal that the Purcell enhancement factor 
remains finite even in the absence of losses, and it depends on the emitter size. 
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I. INTRODUCTION 

Purcell effect is the enhancement of the spontaneous 
emission for a source placed in the resonant cavity as 
compared to that in vacuum jl|]. The engineering of the 
radiative lifetime is now extensively studied in a variety 
of different systems including metallicparticles , mi- 
crocavities p43|, and metamaterials [9l4l3l|. 

The huge Purcell factor for a point dipole embedded 
in the so-called hyperbolic medium has been reported 
in Ref. Ihl . This system, namely an uniaxial medium 
where the transverse e xx = e yy = e± and longitudinal 
e zz = En dielectric constants have the opposite signs, 
is characterized with the hyperbolic isofrequency con- 
tours in the wavevector space [l5l - [l7| , see also the insets 
in Fig. [1] Wave propagation and refraction in the hy- 
perbolic medium reveals its unusual optical properties, 
as compared to the uniaxial medium with the ellisoidal 
isofrequency surface [l5l [l8l - l20j . The radiative rate of the 
point dipole in such a medium diverges, and it remains 
finite only due to the inevitable losses [Tll - ll3l |. 

In this work, we consider a finite-size light source, such 
as a quantum dot, placed in a homogeneous uniaxial 
medium with £m£j_ < 0. We demonstrate that for the 
spatially distributed source the radiative rate does not 
diverge even for vanishing losses, but it depends strongly 
on the source size instead. The maximum enhancement 
of the Purcell factor can be roughly estimated as (A/a) 3 , 
where A is the light wavelength is vacuum, and a is the 
characteristic size of the source. 

Importantly, the hyperbolic medium employed in our 
calculations is not only a hypothetic theoretical model. 
It appears as an effective medium in the homogenization 
of the artificial photonic structures — metamaterials, the 
layered structures consisting of alternating dielectric and 
metallic slabs (2l| - |23| , as well as the structures crated by 
a mesh of metallic wires [24l427l |. In these metamaterials, 
one also have to take into account the stron g sp atial dis- 
persion of the effective dielectric constants [28l l29j and 



the excitation of plasmons [3CJ-[32j. 

The paper is organized as follows. In Sec. [TT] we in- 
troduce our model and outline the calculation technique. 
The results of calculation are summarized in Sec. IIII1 
whereas Sec. IIVI concludes the paper. Details of some 
calculations are given in Appendix [X] 



II. MODEL 

We consider a spherical light source (e.g. a quantum 
dot) embedded into an anisotropic homogeneous medium 
characterized by the dielectric tensor e. Equation for the 
electric field reads 



rot rot E = q„D , 



(1) 



where q a = uj/c, oj is the wave frequency, and c is the 
speed of light in vacuum. The displacement vector in- 
cludes the background contribution and resonant polar- 
ization of the emitter P: 

D = eE + AtiP , (2) 

and the nonzero components of the dielectric tensor are 

&xx ^-yy — ^-L? &zz — £|| • (3) 

We write the phenomenological material equation for the 
polarization as [33|, l34j 



ri 2 $(r) 

h{ujQ — uj) 



J d 3 r'$>(r')E(r') 



(4) 



Here ujq is the resonance frequency, d is the effective ma- 
trix element of the dipole moment of the emitter, and the 
function <&(r) characterizes the spatial distribution of the 
emitter polarization. In what follows, we use $(r) in the 
simple Gaussian form 



%/2exp(-r 2 /2a 2 ) 
47r 3 / 2 a 3 



(5) 
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so that J d 3 r$(r) = 1. Equation Q is similar to the ma- 
terial relation for the excitonic polarization of the semi- 
conductor quantum dot, see Ref. I3"H 

The radiative lifetime t is related to the complex eigen- 
fre quency u> of the homogenous system Eq. (QJ-Eq. (|4]) 
as I34H36II 



— = —2 Imw . 
T 



To find r we apply the Fourier transform, 

d 3 k 



E(r) 



-E k e 



Ikr 



and obtain 

[qleE k -k 2 E k +k(k-E k )} = - 
where 



47r^d 2 $ fc f d 3 k' 
h(uj - u) J (27r) E 



(6) 



(7) 



rE k ,^ k , 



$ fc = / d 3 re" ifer $(r) 



= -k 2 a 2 /2 



(8) 



(9) 



In the derivation of Eq. ©, we took into account that 
the function depends only on the absolute value of 
the vector k. Equation ([5]) can be rewritten as 



E k = -^* k G k A 



(10) 



UIq — UJ 

where we introduced the Green function in the fc-space 

Gk = (M fe ) -1 , M k<a = qoSa/3 - k 2 5 a p + k a kp (11) 

and defined a new variable, A = J E k $ k d 3 k/ (2tt) 3 . Mul- 
tiplying both the parts of Eq. © by and integrating 
over k, we obtain the matrix equation for the complex 
eigenfrequencies w, 



4nq 2 d 2 



d 3 k 



(w - lo q )A = RA, R = J j^^tGk ■ (12) 

We note that the matrix R in the right-hand side of 
Eq. (1121) generally depends on the frequency uj. However, 
we are interested in the weak coupling regime, when the 
interaction of the emitter with light can be treated as a 
perturbation [33], and we set R(uj) = R(ujo) in Eq. (IT2l . 



Taking into account that the matrix R is diagonal due to 
the symmetry of the problem, we obtain the spontaneous 
emission times 



1 _ 8Trd 2 q 2 f d 3 k 



(27T) 



ImG a! fe$^, a = x,y, z. (13) 



The times t x = r y and t z describe the decay of the source, 
initially polarized in the plane xy and along the z axis, 
respectively. To find the decay rates, one should substi- 
tute the explicit expressions for the Green function, 



Gk.zz — 



G k xx 



1 



l-k 2 /(q 2 e ± ) 



% - k ±/ £ \\ 



(14) 
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k 2 



k 2 x [l-k 2 J{q 2 E^)] 
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k 2 /e±_ 
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k 2 /e ± 



*x/ e l|-"-||/ 

into Eq. P^|) . As follows from Eq. (IT41 . the axial dipole 
couples both with TE (ordinary) and TM (extraordinary) 
waves, while the orthogonal dipole couples only with TM 
waves. The mode dispersion can be determined from the 
poles of the Green functions and is illustrated by insets 
of Fig.CQV) and (b). 

The convergence of the integrals (|13|) is assured by the 
rapidly decaying function Thus, the cutoff is natu- 
rally provided by the source size, fc max ~ 1/et, similarly 
as it happens in the nonrelativistic theory of the Lamb 
shift [38 1 . We note, that the realistic metamaterial such 
as the wire medium is characterized by the lattice con- 
stant clq. If the emitter size is smaller than the spacing 
between the wires, a < ao, our approach is not appli- 
cable, and the cutoff is provided at fc max ~ 1/ao- It 
was shown in Ref. l39l . the enhancement of the density of 
TEM modes in the wire medium as compared to the TM 
modes in vacuum, is of the order of l/(<7o a o) 2 - This value 
provides the estimation of the Purcell factor for the wire 
medium. 



III. RESULTS AND DISCUSSIONS 

The integral in Eq. ([TBI can be readily calculated nu- 
merically. In case of the source size smaller than its wave- 
length, one can also obtain explicit analytical expressions 
(see Appendix [Al for more details), 
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(15) 



where 



e ||/ £ J 



Eqs. (fT5j) present the central result of this work. They 
are valid for an arbitrary complex values of £11 and e± 
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FIG. 1: (Color online) Purcell factor relative to vacuum as a 
function of (a) eji for e' ± — 1 and (b) e' ± for eji = 1. Solid 
and dashed lines correspond to the dipole oriented along the 
x and 2 axes, respectively. The insets show schematically the 
isofrequency surfaces in the fc-space for £j_ e '| < 0- Calculation 
was performed for qoa = 0.1 and e'u = e" = 0.1. 



FIG. 2: (Color online) (a) Purcell factor as a function of e'u 
for e" = e" = 0.01 (dashed curve), e'_[ = eff = 0.1 (dotted 
curve), e" = e|f = 0.5 (solid curve) and qoa=l. (b) Same as 
in (a) but for qoa = 0.2 (dashed curve), qoa = 0.1 (dotted 
curve), qoa = 0.05 (solid curve) and e" = e" = 0.1. The 
dipole is oriented along the z axis. 



provided that <Zo a \/k|iT ^ 1> Qo<i\/\s±\ "C 1. The ex- 
perimentally observed decay kinetics of the emitter will 
be determined by the excitation conditions. In case when 
the direction of the dipole moment is fixed and makes the 
angle 9 with the symmetry axis z, the decay will be bi- 
exponential with the initial slope given by 



1 



; e sin^ e 



-(e) 



In the isotropic medium, where e± = e 
(fT5|) reduce to 



(16) 

= k all the rates 



1 

T 



3 j2 



3ft 



Re- 



'3h^/ira'' 



Re- 



(17) 



In the transparent medium (k" = 0) the first term in 
Eq. (|17[) reduces to the textbook result for the sponta- 
neous emission rate (4p| |. The second term describes the 
energy losses due to the heating of the medium [U [i3 
similarly as for a dipole placed in the pore in metal [5|. 
This term controls the decay rate, when the real part of 
the dielectric constant is negative. 

In the anisotropic medium with vanishing losses, i.e. 
e" = e " _> Eqs. JT5]) reduce to 



wi°\ 



(1) 



ir, 



(3) 



ir, 



(i) 



90 a 



(3) 



(9oa) 



s' ± > 0, ej| > 
r , e', > 0,s[ < 



(18) 



q a 



0, 



s' x < 0, ej| > , 
e' ± < 0, e|| < , 



where a = x,y, z, 



3fi 



1 ^(0) = 4 ^ 2 ( £ ll+ 3£, - L ) 
3/iy£ ± 



(19) 



^wr\ q ^(2\s' ± \ + ij 

2^,1 + 1^1)3/3 

v^kjll 2 

4ft^l(l4l + kll) 3/2 ' 



V^lTgo^ 2 ^(3) = _ 

2^(^1 + 1^1)3/2' x 4^^(141 + 1^1)3/2 



and W y = Wjc 



Eq. (11 8[) clearly demonstrates nonanalytical behavior 
with en and ex- When both constants are positive, we 
are dealing with traditional uniaxial dielectric and tran- 
sition rates do not depend on the dipole size. For £|| < 
and e±_ < all radiative rates are zero since the waves 
in such medium are evanescent and do not carry energy 
away from the source. The most interesting regime is 
realized when longitudinal and transverse dielectric con- 
stants are of the opposite sign, £n£j_ < 0. In this case 
the emission rates are governed by the terms oc 1 jo? and 
cx 1/a in (|18p . i.e. are determined by local field effects. 
Counterintuitively, in the regime when £||£j_ < the local 
field is extended in the whole space [TtJ and thus controls 
the radiative emission. The detailed analysis of this pe- 
culiar field pattern will be presented elsewhere. 

The results of numerical calculation of the transition 
rates based on Eqs. (fl"3)) and (jT4j) are summarized in 
Figs. [1] to [31 The rates are normalized to their values 
at £j_ = £|i = 1, which yields the Purcell factor with re- 
spect to vacuum. Figure [1] shows the Purcell factor for 
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FIG. 3: (Color online) Frequency dependence of the Purcell 
factor in the medium with e± = 1 and en given by Eq. (|20|l . 

(a) Functions eji(w) (solid curve) and eji (u) (dotted curve). 

(b) Purcell factor for the dipoles oriented along x (solid curve) 
and along z (dotted curve). Calculation was performed for 

qo a = o.i, e| 0) = i, r/r = 0.03. 



the different dipole orientations as a function of (a) e|, 
and (b) e±. In agreement with Eqs. (fT5)l . the rates dras- 
tically increase when the real part of one of the dielectric 
constants becomes negative. Interestingly, the largest en- 
hancement in Fig. [IJa) is achieved when ej, is negative 
but small, i.e. when e'« < and |ej,| < 1. In agree- 
ment with this result, the leading terms in Eqs. (fT5)l for 
the transition rate are proportional to -y/^jT/ ( a3 l 3 ^ 2 ) 
for |ej| | 3> 1. This justifies additionally the fact that the 
observed enhancement is the local field effect, because 
for large values of |eji| the local field is screened, and so 
the effect is suppressed. Similar analysis applies for the 
dependence of the transition rates on e' ± , see Fig. QJb). 
We notice that the analytical results (fT5j) describe all the 
curves in Fig. Q] with a precision better than 5%. 

Figure [5] shows how the Purcell factor dependence on 
£|| changes (a) with losses and (b) with source size a. 
From Fig. [2^ we conclude that the losses smear the non- 
analytic behaviour of Purcell factor when ej, crosses zero 
and reduce maximum value of Purcell factor. On the 
other hand, in the regime when ej, > 0, e' ± > 0, the losses 
lead to the growth of the decay rate. This is similar to 
the isotropic case, Eq. (|17l) and is related to the heating 
of the medium by the emitted field. Fig. [5p shows, that 
the Purcell factor is very sensitive to the dipole size. It is 
quickly suppressed when the size increases, in agreement 
with 1/a 3 and 1/a terms in Eqs. (fT5)l . 

Figure [3] illustrates the frequency dependence of the 
Purcell factor in the medium where 



£||M 



-(0) 



u>q — uj — ir 



(20) 



and £j_ = 1. Comparing Fig. [3ja) and Fig. [3fb), we 
observe that the largest enhancement is achieved in the 
spectral region where gj, (uj) is negative, but small, which 
agrees with our analysis of Fig. [T] As a result, the posi- 
tions of the maxima of the curves in Fig. [3]Jb) are blue- 
shifted from the resonance energy uiq. 



IV. CONCLUSIONS 

We have developed the theory of the Purcell effect for 
spherical dipole emitters embedded in homogeneous uni- 
axial media, taking into account a finite size of the emit- 
ter and losses in the surrounding medium. We have ob- 
tained analytical expressions for the decay rates in the 
case when the emitter size is much smaller than the wave- 
length of radiation. We have revealed that, when the real 
parts of the longitudinal and transverse dielectric con- 
stants e'| and e'j_ are of the opposite sign (i.e. for the 
so-called hyperbolic media), the radiative decay rate de- 
pends strongly on the emitter size, and it diverges when 
the size vanishes. This enhancement is related to the 
peculiar nature of the local field in such systems, which 
spatially extends to infinity. The largest Purcell factor is 
achieved when < and the absolute values of the 
dielectric constants are much smaller than unity, since 
the screening of the local electric field in this case is min- 
imal. Our theory has been developed for a homogeneous 
medium, but it can also provide a qualitative insight into 
the problem of the spontaneous emission in metamateri- 
als. 
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Appendix A: Analytical expressions for the decay 
rates 



In this Appendix, we present the details of the deriva- 
tion of Eq. (15) . 

First, we substitute Eqs. (TUT) into Eq. (fT5)) and intro- 
duce the spherical coordinates (k,6,(p) in the fc-space. 
Integration over the azimuthal angle dip and over dk can 
be performed analytically, and yields 
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'rad 



, n ( X(0) sin 9 cos 2 9 \ ; 2 (0)sin d ix 5/2 (6>) sin 3 9 
Im I d9\-—= — H 1 g e 



( q oa)\(B) [ X + erf [ igofflV ^)]] | , (Al) 



^ad ^rad ft J I 2 40TI 



x(0) sin 3 



7r£i|e ± (q a) 2 



sinflfeii — -ysin 2 9) iy 3 / 2 (#) sin # , n2 , a \ r * n o 1 

+ -L| i + X , U e-fa°°) *M 1 + erf i go aV^)] [e„ - X sin 2 0] j , A2 



where 



1 sin 2 9 cos 2 9 



(A3) 



and the error function is defined as 
erf(x) = 

We consider the case when the source size is very small, 
so that the condition 



Qoa\y/x(9)\ < h qoa\y/eI\ < 1, goaly^l < 1 

(A4) 

is satisfied for all values of 9. If the real parts of the di- 
electric constants have the same sign, the conditions (|A4|) 



are easily satisfied for small q^a. However, if £m£j_ < the 
quantity x($) may vanish. In this case, the first condition 
(|A4[) will be still satisfied provided the imaginary parts of 
the dielectric constants are sufficiently high. Under the 
conditions (|A4[) the exponential functions in Eqs. (|A1[) 
can be replaced by unity, and the error functions can be 
neglected. 

After that simplifications, the integration over 9 can 
be performed analytically, and it gives Eqs. (|15[) . Our 
numerical analysis shows that Eqs. (| 15[) hold even when 
e" , e'u vanish, provided that the last two conditions (IA4|) 
remain valid. In this case, the terms in Eqs. (IA1|) and 
(|A2|) proportional to exp[— (q a) 2 x(&)] are rapidly oscil- 
lating, so that their contribution to the integrals becomes 
small. 
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